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Abstract. The geometric spreading in plumes and in the interplume region in coronal holes is 
calculated, using analytic and numerical theoretical models, between 1 .0 and 5.0 R ^ . We apply a two- 
scale approximation that permits the rapid local spreading at the base of plumes (//) to be evaluated 
separately from the global spreading (f g ) imposed by coronal hole geometry. We show that // can be 
computed from a potential-field model and f g can be computed from global magnetohydrodynamic 
simulations of coronal structure. The approximations are valid when the plasma beta is small with 
respect to unity and for a plume separation small with respect to a solar radius. 


1. Introduction 

The geometric spreading along a stream or flux tube in a coronal hole is generally 
described in terms of the area, A(r ), or the spreading (or expansion) factor, /(r), 
of a nearly radial, infinitesimal flux tube under steady-flow conditions. In spherical 
coordinates, where Ay is the area of the flux tube at its base, at radius ry, f(r) is 
defined in terms of ^4(?') through the equation 

A(r) = (r/r 0 ) 2 /(r)A 0 . (1) 

Here we show what theoretical predictions are for the variation of f(r) in coronal 
plumes and coronal holes. 

The interest in f(r) derives from its importance for the dynamics of the expand- 
ing solar wind. Kopp and Holzer (1977) introduced the spreading factor formal- 
ism to show that the geometric divergence of coronal holes can impose dramatic 
changes on steady state solutions to the solar wind equations relative to the solu- 
tions in which /(r) = 1. Recently, the same formalism has been used to model 
the flow in coronal plumes, in coronal holes between plumes (the ‘interplume’ 
flow), and to undertake studies of complex physical processes that are generally 
beyond the reach of multidimensional MHD models, such as multifluid flows and 
anisotropic temperatures (e.g., Habbal et aL , 1995; Wang, 1994; Hu, Esser, and 
Habbal, 1997). There have also been recent studies using SOHO instruments that 
place empirical constraints on f(r) (DeForest et al , 1997; Poletto et al , 1997). 
What has been absent are theoretical models of coronal spreading factors. This is 
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despite the existence of MHD models which implicitly contain a partial description 
of the spreading factors, but for which f(r) has not been specifically derived. For 
this reason, we describe here what we believe are typical theoretical predictions 
for spreading in coronal plumes and coronal holes between 1.0 and 5.0 Rq. The 
immediate purpose is to show what physical models, based on our present under- 
standing of conditions in coronal holes, predict as an appropriate choice for the 
/(r ) in one-dimensional (ID) flow models. 

There are two parts to the calculation of f(r). The first is the spreading at the 
base of plumes, up to heights comparable to the typical plume separation distance, 
D, where D k 35 000 km (Suess, 1982; DeForest et al, 1997). This calculation 
accounts for the location of plumes over concentrations of magnetic flux relative 
to the background distribution of flux in the interplume region (DeForest et al, 
1997). The flux from these concentrations spreads out with increasing height, as in 
the spreading above the chromospheric network described by Gabriel (1976), until 
the field is approximately uniform. Because the magnetic field is relatively strong, 
this happens by a height comparable to the distance between flux concentrations 
(Newkirk and Harvey, 1 968; Del Zanna, Hood, and Longbottom, 1 997). The second 
part of the calculation is the spreading imposed by the geometry of the coronal 
holes in which plumes lie. This occurs over a length scale that is large relative to 
the typical distance between plumes. Therefore, plume spreading is most simply 
thought of, and can be mathematically approximated as a ‘two-scale’ problem in 
which the local spreading at the base of plumes is computed separately from the 
larger scale spreading defined by the geometry of the coronal hole. This approach 
can be made formal by assuming the separation can be made and then testing it 
a posteriori, just as done in, for example, boundary layer two scale approximations 
(Hinch, 1991). 

Before calculating /(r), we will describe two general, and basically intuitive 
results for how the spreading of the plume and interplume regions behave. These 
derive directly from the observation that the plasma beta, ft = (\6wnkT/ B 2 ), is 
much smaller than unity, ft <C 1.0, throughout coronal holes, up to heights of at 
least 5 7?y. The first result is that the spreading near the base can be computed 
from a potential field model. This has already been shown by Newkirk and Harvey 
(1968), Suess (1982), and Del Zanna, Hood, and Longbottom (1997) and needs 
little or no additional proof. The second is that higher up in coronal holes the 
geometric spreading factors along streamlines inside plumes and in the adjacent 
interplume region vary together, so that if the plume area doubles between 2 and 
5 7?0, then so must the interplume area. This follows directly from ft <C 1.0, but 
requires additional discussion since the field there is not a potential field. 

To demonstrate the important point that ft < 1 in coronal holes, consider 
conditions at the base of plumes, at the top of the transition region at a height 
of ~ 7000 km above the photosphere. There, B = 0[ 1 0-1 00] G (DeForest 
et al, 1997), T = O[10 6 ] K, and n 0 = 0[10 7 -10 9 ] cm" 3 (Ahmad and With- 
broe, 1977). Therefore, ft = 7 x 1 0^ 2 -7 x I O' 6 and is easily «; 1. Outside 
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plumes, in a coronal hole, both the magnetic field and the density are slightly less, 
D = O[5-10] G (DeForest etal., 1997), and n = O[10 7 - 10 8 ] cm" 3 (Ahmad and 
Withbroe, 1977), so that/3 = 3 x 10"" 2 -7 x 10~ 4 . Conditions higher up in coronal 
holes can be estimated by extrapolating inward from solar wind observations. Suess 
et al. (1996) estimate that ft ss 3.0 at 1 AU from Ulysses measurements. Using the 
constant velocity approximation and taking the temperature to be T = O[I0 6 ] K 
gives /3 < O[10~ 2 ] at 5 Rq. Even if the temperature were as high as 0[ 1 0 7 ] K, as 
reported recently for the proton kinetic temperature by Kohl, Strachan, and Gardner 
(1996) and Strachan et al. (1997), it would still only give (3 < O[10~'] at 5 /?©. 
Because fl decreases with decreasing height, it can safely be assumed that /3 <§C 1 .0 
below 5.0 i?0. 

The two-scale problem and associated definitions are described in Section 2, 
the rapid spreading near the base of plumes will be calculated in Section 3, and the 
spreading in coronal holes and in plumes above the rapid spreading region will be 
calculated in Section 4, followed by a brief discussion and summary of the results. 


2, Two-Scale Analysis 


Given the definition of the spreading factor in ( I ), the magnetic field variation along 
a nearly radial field line can be written as 


B{r) = Bo 



( 2 ) 


In coronal holes f(r) is often assumed to increase from unity to some number of 
O[10] by 5 Rq, and change slowly, if at all, after that (e.g., Kopp and Holzer, 1977; 
Hu, Esser, and Habbal, 1997)). As described in the Introduction, it is assumed 
here that there is superimposed on the global coronal hole spreading analyzed in 
those earlier studies a local spreading at the base of plumes which changes rapidly 
between the bottom and some height 5R^ that is of the same order as D, the plume 
separation distance. With this assumption, we can mathematically separate the two 
scales by making substitution 


f{r) = fi{r)f g {r) ■ (3) 

fi(r) is the local spreading below d/?©, varies rapidly at these heights, and is 
constant above 35 000 km. f g (r) is the global spreading, which varies much more 
slowly than //(r) below 35 000 km. 

Using (3), Equation (2) becomes 

RcA 2 1 


B(r) = Bo 


r 


fl{ r )fgi r ) 


(4) 
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The solenoidal condition, V • B, then separates into two equations according to the 
two scales to give 


» d(r 2 /(r)) d/< d/i 

r- dr 9 dr 9 dz ! 


R& < e < 7?.0 + SR^ , 


I d(r 2 f(r)) = 1 d(r 2 / g ) 

r 2 dr 1 r 2 dr 


< r < oo , 


(5a) 


(5b) 


where 2 is a local Cartesian coordinate defined such that e*||e r and z — 0 at 
r = 7?0. Equation (5a) has reduced to the Cartesian form because the thickness of 
the layer is small compared to its radius of curvature. Therefore the geometry in 
this thin layer is Cartesian if SRq <§C R^ and the magnetic field variation across 
this layer is given by 


B(z) = Bo/ft(z) , 


(6a) 


while globally the field variation along a nearly radial magentic field line is given 
by 


B{r) = B(z = SRu) 


Rq + SR^V 


1 

f<A r ) ' 


(6b) 


Now consider how well the magnetic field is simulated by a potential field for 
0 < z < dRy. This can be evaluated by analyzing the momentum equation as 
written in dimensionless form (Del Zanna, Hood, and Longbottom, 1997): 


(V X B) X B = (0/2 )[M 0 2 p(V ■ V)V + Vp + g 0 pr- 2 e r ] , (7) 


where the quantities have been nondimensionalized against typical coronal base 
values, (3o ; T 0 ,rc 0 ), Mo — Vq/Ko (Ko the reference sound speed), and (jo ~ 
SR^GM^rnp/ilkToR^). The V operators have been scaled by 5R^ since the 
independent variables may change by as much as an Of 1 ] amount across the thin 
layer at the base of plumes, while the radius in the gravity term is correctly scaled by 
R^j. Consequently, (6) takes on a slightly different form than given by Del Zanna, 
Hood, and Longbottom in that the definition of involves the ratio SR^/R^ . With 
this scaling, all of the dimensionless independent variables and operators inside the 
square brackets in (6) are assured of being 0[1] or smaller. 

Equation (7) now provides an easy way of estimating the relative importance 
of the various terms by estimating the values of /?, Mo, and go. We have already 
shown that (3 <C 1 .0. Since the flow is subsonic below 2-3 7?^, then Mq 1.0 as 
well. Finally, go — Ofl], so each complete term inside the square brackets in (6), 
and not just the individual variables and operators, is of Ofl] or smaller, showing 
that the left-hand side of (6) is of O [ft]. So, to Ofl], (V x B) x B — 0. Since 
we impose no currents flowing through the lower boundary, the magnetic field 
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Figure /. Distribution of flux concentrations on the plane at the base of plumes in the ‘mesa’ model 
(see text). 


is therefore a potential field to within a formal precision 0[(3], This is the point 
made by Newkirk and Harvey (1968) and Suess (1982). Moreover, Del Zanna, 
Hood, and Longbottom (1997) demonstrated the validity of the approximation in 
a quantitative way by showing the small size of first order effects of the presence 
of a plasma. 


3, Spreading at the Base of Plumes 

For the present purpose, the important property of the magnetic field at ^ = 0 is that 
there are quasi-regularly distributed unipolar flux concentrations and a background 
magnetic field, as most recently found to be the case using MDI observations 
(DeForest et al , 1 997). The example we use has the flux concentrations distributed 
on vertices of a rectangular grid as given by the expression 

Bo{x, y) = g(x, y)( 1 - co) + q> , (8) 

where cq is a constant giving the ratio of the background average field strength to 
the peak field strength in the flux concentrations. The function g(x , y) is arbitrary, 
chosen in such a way as to simulate the shape, areal distribution, and density of 
magnetic flux concentrations on the z = 0 plane. The magnitudes of the field 
strength and g(x, y) are normalized to unity. An example of Bo(x, y) is shown in 
Figure 1. 

The magnetic field above z — 0 can then be shown, using simple potential field 
theory, to be given by B = — V$, with 

oo 

$ — -CQZ + b m b n exp[— 7 mi „z] cos(m7r:c) cas(nwy) , 

m, 7i 


(9a) 
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X 


Figure 2. The total magnetic field strength at the base of plumes (z = 0) for the ‘mesa’ and ‘spike 
distributions. 


7m, n = 7r(m 2 + n 2 )' /2 . (9b) 

This is just the three-dimensional (3D) analog of the 2D calculation described by 
Suess (1982) (except for the correction in (9b) of a factor of 7r error that was noted 
by Del Zanna, Hood, and Longbottom (1997)). The coefficients in (9a) are derived 
by doing a Fourier decomposition of the flux concentrations using the expression 


2 2 


d m ,n = J j a(x, V) cos(m7r:r) cos(mTry) dx d y , 


(10a) 


o o 


with (for (ra,n) 7^ (0,0)) 

b m b n = (1 ~ C ° )fim,n ■ (10b) 

7m,n 

Because this is a potential-field model, the only length scale is the dimensionless 
ratio of the full width at half-maximum (FWHM) of the flux concentration to the 
separation between flux concentrations - there is no dimensional scale. Similarly, 
since the equations are linear in field strength, there is no natural field strength 
and the only amplitude scale is the ratio co. Dimensional values will be introduced 
later. The field exponentially approaches a smooth value with increasing height, 
reaching the average value g( 1 — Co) + Co for z 2> the FWHM of the concentration 
(<7 = the average of g(x, y ) over the ( x , y ) plane). The asymptotic value of the 
spreading factor can be determined by the percentage of the total flux in the flux 
concentration and does not require specific evaluation of f{r). 

The new calculation we do here is to show the detailed profile of fi{z), which 
we do for two cases with different FWHM. These two cases are shown in Figure 2, 
where the total field strength (Equation (8)) is plotted at the base of the plume along 
the ;r-axis (the y-axis is identical). The broader, flatter distribution has a FWHM of 
~ 0. 14 and the thinner, more peaked distribution has a FWHM of ~ 0.08. In both 
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(a) (b) 




Figure 3 . (a) The total field strength at the base (z = 0) and at z = 0.5 for the ‘mesa’ model of a 
plume magnetic flux concentration. The peak to background ratio is 20 at the base but is only 1.05 
(5%) at z = 0.5. (b) Field lines in the x - z (or y - z) plane for the ‘mesa’ model of a plume. 


cases, Co = 0.05, meaning a 5% background field. We will call these distributions 
the ‘mesa’ and ‘spike’ distributions, respectively, and use them to illustrate the 
effects of varying shape and FWHM on the solution. The mesa was illustrated in 
Figure 1 . 

Equations (9) and (10) were solved for the magnetic field above jz = 0 and 
Figure 3(a) shows the resulting field variation at z = 0.5 and z = 0.0 for the 
mesa model. It is seen that the field is essentially smooth by z = 0.5. As noted in 
the Introduction, plumes and flux concentrations are typically 35 000 km apart so 
z = 0.5 corresponds to a height of 17 500 km. This meets one of the criteria for the 
two scale analysis - that the field be smooth by a height 5R® <C R& where we will 
henceforth conservatively take 5R& = D = 35 000 km. Figure 3(b) graphically 
shows the shape of the field lines in the x-z plane, illustrating the typical spreading 
at the base and the rapid approach to nearly vertical field lines by z = 0.5 that is 
implied by the result in Figure 3(a). Note that the field lines are already rapidly 
diverging at z = 0. 

The spreading factors are found from Equation (6a) by computing the variation 
in field strength along a given field line and then dividing that strength by the height. 
Values are plotted in Figure 4(a) along the field lines labeled A-F in Figure 3(b). 
The spreading factor has a large positive gradient at z — 0, increasing quickly to 
its maximum value, f° r field lines A-C. Field line C starts at x ~ 0.15, 

well down on the flank of the concentration (see Figure 3(a)), but still reaches 
approximately the same maximum value of fi %ma x reached along field lines A and 
B. Field lines E and F are outside the flux concentration and exhibit values of fi{z) 
slightly less than unity, as must be the case for magnetic flux conservation. 

The variation of fi(z) is seen from Figure 4(a) to depend only weakly on 
position within the flux concentration. Therefore, we can show the characteristic 
dependence of fi {z) on FWHM and co by plotting just the variation along field line 
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Figure 5. The magnetic field lines in a quasi-steady model for the global MHD structure of the solar 
corona with a dipolar magnetic field, a momentum and heat sources, and thermal conduction (W97). 
The axis of symmetry is horizontal and the global topology is similar to that for any model with a 
heat source. Models without the heat source (SSW, W95) do not have the sharp cusp at the top of the 
streamer. 


relatively slow and dense solar wind which poorly simulates actual coronal hole 
conditions. W95 is similar to SSW except that (3 = 1.0 and it uses a different 
numerical technique. The density is even higher in the open fields at the poles than 
in SSW. W97, using the numerical technique of W95, adds heating and momentum 
sources and thermal conduction. Furthermore, the boundary conditions depend 
on polar angle. The consequence is that W97 simulates empirical coronal hole 
densities, temperatures, and flow speeds in the open regions at the poles while still 
keeping j3 rs 0.5 at the equator, beneath the streamers. The magnetic field lines for 
W97 are shown in Figure 5, where the field line geometry is similar to that for all 
models with volumetric heating ( Steinolfson, 1 988; Suess et al., 1 996; Wang et al., 
1997). W95 and SSW, without volumetric heating, have more rounded streamers 
than those shown in Figure 5. The difference is due to the slow expansion of the 
streamers and leakage of slow solar wind out of the tops of the streamers when 
there is volumetric heating, a process which is called streamer evaporation (Suess 
et al., 1996). 

The spreading factor, f 9 (r), was computed using Equation (6b) along two field 
lines for each example; one near the center of the open field region (9 = 0°) and one 
near the boundary of the streamer. First the coordinates along a magnetic field line 
were computed from the base of the corona upward from the numerical solution, 
then the field strength along the line was inserted into (6b) to give f 9 {r). Because 
the streamer width varies between the models, the starting polar angles for the field 
line near the boundary of the streamer vary. The spreading factors between 1 .0 and 
5.0 are shown in Figure 6. 
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Figure 6. Spreading factors from the MHD global coronal models SSW (a), W95 (b), and W97 (c). 
These models are described in the text. In each panel, the spreading factor along a field line near the 
center of the coronal hole (0 < 2.25°) is shown as a solid line and that along a field line near the 
edge of the coronal hole (0 > 38.25°) is shown as a dashed line. The indicated angle for each field 
line is the angle at the footpoint, at 1.0 = R&. 


Figure 6 can be summarized by noting that although there are broad similarities 
between the three models, there are several specific differences. The similarities are 
that the spreading varies smoothly with radius and /^ max < 6.0 in all the cases - 
although the W95/9 = 1° result may become larger than 6.0 at larger radii. Potential 
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Figure 8. Comparison of the global spreading factor in the B = 1° direction from Steinolfson, 
Suess, and Wu (1982), the plume spreading factor from Wang (1994), and the Habbal et ah (1995) 
interplume spreading factor between 1.0 and 5.0 R&. 


using the exponential variation suggested by Equation (9). This would apply to 
the spreading in plumes. A corresponding, but much smaller, decrease in //(r) to 
values less than unity should occur in the interplume region in this thin layer at the 
base of plumes. 

A more serious conflict is found in the choice by Habbal et al (1995) of 
different rates of change for f(r) in the plume and interplume above heights of 
35 000 km. The physical arguments and theoretical models described here strongly 
suggest that the relative changes in the spreading factors should be the same. 
The choice by Habbal et al (1995) was based on a reported empirical result 
that the spreading factor in plumes w^as near unity, based on the observed widths 
in brightness of plumes using SPARTAN 201 (Guhathakurta and Fisher, 1995). 
The only explanation for this discrepancy is that the widths in brightness do not 
necessarily correspond to the width that would be computed by following a specific 
stream or field line, since the two are not necessarily the same (Suess, 1982). 

More generally, the spreading factors utilized in these two ID models are well 
within the range of possibilities for f g (r) derived from global MHD models - 
the SSW curve for f g (r) is overlaid onto Figure 8 to illustrate this point. This is 
reassuring since the numbers used in the ID models were derived from estimates 
based on the empirical geometry of boundaries of coronal holes (Munro and Jack- 
son, 1977) which, by the arguments presented here, should be typical of the values 
throughout coronal holes. It means that the values for fg(r) derived from all the 
MHD models are consistent with values derived empirically for large polar coronal 
holes near sunspot minimum. 
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6. Summary 

The geometric spreading in coronal plumes and coronal holes has been derived 
here from theoretical models, based on the observation that /I « 1.0 throughout 
coronal holes. The calculations were done in two parts. First the spreading below 
SRq — D was found, where D = 35 000 km is the typical distance between 
plumes or magnetic flux concentrations lying at the base of plumes (there are 
probably many more flux concentrations than plumes). This spreading was found 
using a potential field model. It showed that plume spreading is ~ 15 for a 5 % 
background interplume field and that there is a small corresponding convergence in 
the interplume region at these heights. The field is then locally smooth above SRq. 
Second, it was shown that the global spreading above SR^ could be calculated 
using global MHD models and this was done for three different published models. 
The spreading was determined to be similar to empirical spreading found from 
the geometry of coronal holes boundaries, f g ^ m ax ~ 1.5-6 at 5.0 R&. Combining 
these two results produces an overall spreading in plumes of ~ 40 between 1 .0 and 
5.0 R with most of this increase occurring below 5R@. 

An important general result is that plume and interplume relative spreading is 
the same to within a factor of 0[/5] above <LR©= 35 000 km. Plume and interplume 
spreading differs strongly between plume base and <*>/?©, over which distance the 
initially nonuniform magnetic field becomes smooth on transverse scales small 
compared to a solar radius. However, if a plume area doubles between, (e.g., 1.5 
and 5.0 i?©, then the nearby interplume area will also double over this height 
(maintaining, of course, any relative differences between the plume and interplume 
spreading which existed at the height 5R&). This point is made here specifically 
because there are reported differences in empirical relative spreading between 
plumes and interplume in this part of the corona. It will be important to understand 
the source of these differences since the physical premise for the conclusions 
reported here seems strong. 
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